On the ergoregion instability in rotating gravastars 
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The ergoregion instability is known to affect very compact objects that rotate very rapidly and do not possess 
a horizon. We present here a detailed analysis on the relevance of the ergoregion instability for the viability 
of gravastars. Expanding on some recent results, we show that not all rotating gravastars are unstable. Rather, 
stable models can be constructed also with J/M 2 ~ 1, where J and M are the angular momentum and mass of 
the gravastar, respectively. The genesis of gravastars is still highly speculative and fundamentally unclear if not 
dubious. Yet, their existence cannot be ruled out by invoking the ergoregion instability. For the same reason, not 
all ultra-compact astrophysical objects rotating with J/M 2 < 1 are to be considered necessarily black holes. 

PACS numbers: 04.40.Dg, 04.30.Nk, 04.25.Nx 



I. INTRODUCTION 

Gravastars have been recently presented by Mazur and Mot- 
tola 0]] as a new exact solution to the Einstein equations. In 
the original suggestion these are compact, spherically sym- 
metric and non-singular objects, that can be taken to be almost 
as compact as black holes. In this "three-layer" model, the 
high compactness is supported by a de-Sitter core, surrounded 
by a shell of matter, and the exterior vacuum spacetime is, of 
course, that of the Schwarzschild solution. 

According to the original picture proposed in ref. a 
massive collapsing star could go through a phase transition 
when its radius approaches R = 2AI, forming a gravastar in- 
stead of a black hole. Although the dynamical processes that 
would lead to the formation of a gravastar are far from being 
understood and could probably be described as rather exotic, 
the final state can be described by an exact (and fairly simple) 
solution of the Einstein equations. This is what stimulates our 
interest in this matter. 

The original gravastar model has inspired many subse- 
quent works, but always describing stationary solutions (ex- 
cept from ref. 01, as we will see below). A related simpli- 
fied model with an infinitesimally thin shell was proposed in 
ref. JH] and later also generalized in ref. 0]. Several possibil- 
ities for the interior solution have been considered in refs. fH, 
|0] and |0], among others. More recently, a solution for 
electrically charged gravastar configurations was proposed in 
ref. il] and limits on the existence of gravastars from astro- 
nomical data were considered in ref. @|. 

Also an effort has been made to assess the properties of 
gravastars when these are perturbed 03 [H. The result of this 
analysis has lead, among other things, to the evidence that per- 
turbations in spherical gravastars can be used to discriminate 
them from black holes, thus removing one of the most seri- 
ous consequences of the existence of such ultra-compact ob- 
jects Hill . Note that hereafter we define as ultra-compact any 
stellar object whose compactness fi = M/R is much larger 
than that of typical neutron stars of comparable mass and an- 
gular momentum, and which is \i ~ 0.15 — 0.2. This definition 
is inevitably weak and ambiguous, but it aims at focussing on 



the large compactness of gravastars as the most relevant prop- 
erty. Clearly, black holes are always more compact than any 
possible gravastar model (although their compactness is only 
infinitesimally smaller for nonrotating models) but the rele- 
vant point to bear in mind is that gravastars have compactness 
much larger than that of standard stars and comparable (al- 
though smaller) to that of black holes. 

Further expanding on the perturbative analysis carried out 
in ref. fiUl . Cardoso et al. [2] have recently considered the 
properties of perturbed and rotating gravastars and assessed, 
within the slow-rotation approximation, their stability against 
the "ergoregion instability" . We recall that such instability 
affects rapidly rotating and very compact stellar objects which 
do have an ergoregion but do not have an event horizon lfl2l 

a 

The Kerr black hole is a good example of an object that has 
an ergoregion. But an ergoregion can also develop in compact 
stars that are sufficiently rapidly rotating. In this region, the 
relativistic frame dragging is so strong that no stationary orbits 
are allowed. All trajectories of particles in this region must 
rotate in the same direction of the rotation of the star. 

Because of this effect, some particles in the ergoregion can 
be measured by an observer at rest at infinity as having neg- 
ative energy. This happens because this observer measures 
the energy of the particles by projecting their four-momentum 
vectors onto his four-velocity. As the observer is stationary 
(and no stationary trajectories are allowed inside the ergore- 
gion), his four-velocity is outside of the light cone of the par- 
ticles in the ergoregion. This causes some of these particles 
to have their energy measured as negative by the observer at 
infinity. 

The instability occurs then in the following way. From an 
initially small perturbation with negative energy trapped in 
the ergoregion, one can extract positive energy (that leaves 
the star and goes to infinity) by increasing the negative en- 
ergy inside the ergoregion (thus conserving the total energy). 
As the negative energy trapped in the ergoregion increases, 
the star radiates even more positive energy to infinity and this 
process leads to the instability. This process is very general, 
and scalar, electromagnetic and gravitational waves become 
unstable in a star with an ergoregion. 



2 



This paper is dedicated to reconsider the analysis carried 
out by Cardoso et al. [2] and to extend it to a larger space of 
possible models, taking into account the limits on the thick- 
ness of the matter shell and on its compactness. When doing 
this, we confirm the results of Cardoso et al. 12J for their mod- 
els, but also show that the conclusions drawn were excessively 
restrictive. In particular we show that not all rotating gravas- 
tars are unstable to the ergoregion instability. Rather, we find 
that models of rotating gravastars without an ergoregion (and 
therefore stable) can be constructed also for extreme rotation 
rates, namely for models with J/M 2 > 1, where J and M are 
the gravastar's angular momentum and mass, respectively. 

The paper is organized as follows: in Section|II]we briefly 
review our gravastar model and the slow-rotation approxima- 
tion. In Section [III] we develop the equations for scalar per- 
turbations and present the WKB approximation used. In Sec- 
tion [TV] we present our results and analyze the behavior of 
the instability in the space of parameters and in SectionlVlwe 
present our concluding remarks. We use c = G = 1 through- 
out the paper. 



H. A ROTATING GRAVASTAR MODEL 

We start our analysis by adding a uniform rotation to the 
simple fluid grav astar model with anisotropic pressure pre- 
sented in ref. ifllTl - We recall that the use of anisotropic pres- 
sures was introduced by Cattoen et al. Il4ll to remove in part 
the complications produced by the infinitesimal shells in the 
original gravastar model of Mazur and Mottola jj]]. In this 
way, the anisotropic pressure replaces the surface tension in- 
troduced by the matching of the metric in the infinitesimally 
thin shells. Although the use of an anisotropic pressure is es- 
sentially arbitrary, as arbitrary are the equations of state used 
to describe such a pressure, it has the appealing property of be- 
ing continuous and thus of allowing one to build equilibrium 
models without the presence of infinitesimally thin shells and 
thus look more seriously into the issue of stability. 

We write therefore the line element for a finite-thickness 
rotating gravastar with anisotropic pressures within the slow- 
rotation approximation, (namely at first-order in the angular 
velocity in the small parameter 0,/Q.k, where fljc is the Kep- 
lerian limit) as 



ds 2 



z v ^dt 2 



e x ^dr 2 



2 • 2 

-r sin 



r 2 d6 2 ^ 
-uo{r)dtf 



(1) 



together with the energy momentum tensor given by 

T» v = (p + Pt)u"tt* + PtS** v + (p, - ptK*„ , (2) 

where p is the energy density of the gravastar, p r and p t are 
the radial and tangential pressures, respectively. The vector 
it M is the fluid four-velocity, 

u^u^ = -1 , u r = u e = , = nu* , 

«* = [- {att + 2^3*0 + n 2 g ^)] ~ 1/2 , 



and the vector s M is orthogonal to the fluid four-velocity, 

u% = (s t = S<t> = 0) . (3) 



s% = 1 



The functions v{r) and A(r) in the line element (Q]) are 
given in terms of the mass m(r) and the radial pressure p r (r), 
as in the non-rotating case, 



e -x = 1 _ 2w 0) 



2m(r) + 87rr 3 p r 
r(r — 2m(r)) 



(4) 
(5) 



where a prime denotes a (total) radial derivative. As usual, 
we define the gravitational mass within a radius r as m(r) = 
Anr 2 pdr and 



p(r) 



Po , 



< r < n 



ar + br + cr + d, r\ < r < r-2 



(6) 



, r 2 < r 

with the coefficients a, b, c, d given by 

2po 

, _ 3p (r 2 +n) 
p 

_ 6p rir 2 

C - S 3 ' 

p {4 - 3rirf) 



d = 



5 3 



where S = 

Po 



r\ is the "thickness" of the gravastar and 
15M 

' 27r(ri + r 2 )(2r 2 + r x r 2 + 2rf) ' 



(7) 
(8) 
(9) 
(10) 

i 

(11) 



The equation of state used for p r (p) serves here only as a 
closure relation and is therefore chosen to be in the simplest 
possible form, namely a polynomial of the type lflol[T5ll 



Pr{p) 



a _( Q + l) IJL 
Pa 



Po/ 



(12) 



where a = 2.2135 is determined by demanding that the max- 
imum sound speed c 2 at which d 2 p r /dp 2 = coincides with 
the speed of light to rule out a superluminal behavior. Fi- 
nally, the tangential pressure pt is given by the anisotropic 
TOV equation 



r 1 

Pt =Pr + »Pr + + 



m(r) + Airr 3 p r 



r(l — 2m(r)/r) 



(13) 



The function u(r) can be shown to be of first order in the 
angular velocity 51 and it describes the dragging of inertial 
frames. Formally, it can be obtained from the * ^ component 
of the field equations lfl6ll 



Rj = 8tt7\ 



(14) 
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which gives, for our anisotropic case, 



47rr 2 ( i o + p r ) 
r — 2m 



167rr(p +p t ) 
r — 2m 



where 



vj{r) = Q — lu(t) 



w, (15) 



(16) 



for the region of the spacetime interior to the gravastar. It is 
easy to see that eq. (IT3T > is equivalent to eq. (2.25) of ref. iQ]. 
In the exterior region, on the other hand, we have 



w(r) = 



2J 



2J 



(17) 



The definition of w in eq. (fTBT l deserves some attention and it 
corresponds to the difference between il, the angular velocity 
of the gravastar (as seen by an observer at rest), and uj(r), 
which gives the angular velocity of a zero angular momentum 
observer (ZAMO). Therefore, w(r) is the angular velocity of 
the gravastar as seen by the ZAMO. 

As a general result in the asymptotically flat limit, the drag- 
ging must go to zero as co = 2 J/r 3 + 0(r~ 4 ) for r — > oo, thus 
defining J as the angular momentum of the spacetime lfl7ll . 
Demanding now that the exterior equations (fTTT i to be consis- 
tent with this asymptotic limit implies immediately that the 
integration constant J in (T7\ should be identified as the total 
angular momentum of the gravastar. 

A representative behavior of the frame-dragging function 
ui(r) is shown in Fig. Q] for a gravastar with thickness of 
the shell 8/M = 0.4 and compactness /.t = A//V2 = 0.45. 
The frame-dragging is obtained by numerically integrating 
eq. (fIBl . with initial conditions w'(r = 0) = and w(r = 0) 
finite. As shown in the figure, the solution for ui(r) in the 
interior of the gravastar is constant [cf. eqs. (15[ and (fToT ll 
and monotonically decreases outwards. The two integration 
constants SI and J can then be determined by matching the in- 
terior and exterior solutions at the boundary r = r2 ■ Finally, 
in the exterior of the gravastar, ui(r) goes asymptotically to 
zero [cf. eq. (flTt l. 



III. SCALAR PERTURBATIONS AND THE ERGOREGION 
INSTABILITY 

We recall that if a compact relativistic star rotates suffi- 
ciently fast, it will posses an ergoregion, namely a region in 
which the frame-dragging is so intense that stationary orbits 
are no longer possible. Within this region even massless par- 
ticles, such as photons, will be involved in the "drag" and all 
trajectories will rotate in the prograde direction. Mathemati- 
cally, the boundary of the ergoregion is defined as where the 
covariant tt-component of the metric has a zero, i.e. 



9tt 



r 2 L0 2 sin 2 1 



= 0, 



(18) 



and it should be noted that the ergoregion does not need to be 
restricted to the exterior of the compact star and, rather, it can 
also involve regions interior to the star. 




FIG. 1: Typical example for the dragging of inertial frames uj(t), for 
a gravastar with = 0.45, 8/M = 0.4 and J/M 2 = 1 (fi/fijf = 
0.82). We can see that lu — const, in the interior (r < ri) and 
u> — > in the exterior region (r > r-2). 



The main interest in the existence of an ergoregion in a 
compact star stems from the fact that these regions lead to a 
secular instability by means of which any initially small per- 
turbation will grow exponentially in time EL [l2ll . Clearly, one 
expects that nonlinear effects will intervene to limit the growth 
of the instability once this has reached a sufficiently high (sat- 
uration) amplitude which, however, cannot be determined on 
the basis of a linear perturbative investigation. While deter- 
mining such an amplitude is of great physical interest, we will 
here limit ourselves to a simple linear analysis of the problem. 
On the other hand, we will determine not only if an instabil- 
ity can or cannot develop but also, and more importantly, the 
characteristic timescale for the growth of the instability. In- 
deed, as shown already many years ago in ref. 11211 . for ultra- 
compact uniform-density stars, the timescale for the growth of 
scalar perturbations via the ergoregion instability can some- 
times be several orders of magnitude larger than the age of the 
universe. When this is case, the instability grows so slowly 
that the stars are effectively stable even if the linear analysis 
reveals that they are mathematically unstable. 

Hereafter we will concentrate on scalar perturbations and 
within the WKB approximation. There are many references 
in the literature where this method is described, but we sug- 
gest for instance ref. Il8ll . The WKB method (from Wentzel, 
Kramers, Brillouin) is used to approximate the solution of dif- 
ferential equations of the general form 



y"(x) +iv 2 r(x)y(x) = . 



(19) 



where ui 2 ^ 1. As we will see below, the scalar perturba- 
tions that we are considering are described by eq. d22l) . which 
has the same form as eq. ( Tl~9b above. The same equation is 
often found to arise in quantum mechanics, but also in clas- 
sical physics problems. The key point of the approximation 
is the assumption lu 2 ^> 1 (high frequency approximation in 
quantum mechanics). It allows us to consider that, for every 
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u> n 1, then r(x) is approximately constant between two ze- 
ros of y n (x). It is easy to see that this approximation greatly 
simplifies the integration of eq. ( fT9l l. 

Of course, studying the response of a rotating gravastar to 
electromagnetic or gravitational perturbations would be astro- 
physically more interesting and realistic. However, the use of 
scalar perturbations has the important advantage that in this 
case one can decompose the perturbations in spherical har- 
monics and reduce the perturbation equation to a single ordi- 
nary differential equation, which can be integrated with very 
modest computational costs and high accuracy. In addition, 
because the order of magnitude for the growth of the pertur- 
bations is expected to be the same for scalar, electromagnetic 
and gravitational perturbations [this is indeed the case for the 
decay of perturbations of different spins in a Schwarzschild 
spacetime (see ref. lfl9tl for some numerical results)], estimat- 
ing the growth time of scalar perturbations allows for a simple 
and direct extension also to other types of perturbations. 

Bearing this in mind, we next proceed to the study of the 
massless scalar wave equation in the slowly rotating gravas- 
tar background that we introduced in Section [II] We briefly 
review here the basic steps usually followed in this type of 
analysis; such steps can be found in many papers and are re- 
produced here only for completeness. The wave equation 



1 



d 



dx v 



0. 



(20) 



can be separated by using the ansatz 



ij){t,r,6,4>) = xMexp i~ 



x e lat Ye. 



2 v ' A '\ J 1 
— I \dr > x 

r 2 2j J 

(21) 



where the eigenfrequency a is in general complex, with a = 
ay — i jr. As a result, r > leads to an exponential growth, 
while t < leads to a decay. 

In the high m limit (and taking I = m), eq. ( f20l > can be 
written using eq. d2Tb and keeping only the dominant terms 
(i.e. first order in co and 1/m) as 



X. 



m 2 T(r,£)x = : 



(22) 



where m is here the order of the Y( m (6,<f) spherical har- 
monics (not to be confused with the mass function m(r)), 
S = a/m is the negative of the pattern speed of the pertur- 
bation and T is given in terms of the two rotationally split 
"effective potentials" V+ and V-, 



r£ — V 



V± = -u± 



-(E-V + )(£-V_), 

e u/2 



r 



(23) 
(24) 



It is easy to see that the boundary of the ergoregion given in 
eq. ([181 1 coincides, in the equatorial plane, with V+ = 0, so 
that the ergoregion is effectively contained in the region where 
V+ < 0. 

Our stability analysis translates therefore to finding the 
complex eigenfrequencies a of the scalar wave modes that 




FIG. 2: Typical example for the potentials V±, for a gravastar with 
H = 0.45, 5/M = 0.4 and J/M 2 = 1 (Q/Q K = 0.82). The first 
three unstable modes with negative energy "trapped" in the potential 
well are depicted, as well as the points r a , rb and r c for the i = m = 
1 mode. 



satisfy eq. d22b . However, because these frequencies gener- 
ally have the real part being much larger than the imaginary 
one, i.e. , 0?(cr) 3> 9(c) (or, equivalently, ay ^> 1/r), it 
is a reasonable approximation to consider a as essentially a 
real number. It is also more convenient to the problem to use 
S = a/m instead of a itself (see the Appendix for more de- 
tails). 

Furthermore, when considered within the WKB approxi- 
mation, the potentials V± (r) determine four different regions 
which are reminiscent of those appearing for waves trapped in 
the potential well of V+. To illustrate this in more detail we 
show in Fig. |2]an example of the two potentials for some typ- 
ical parameters. In addition, the figure shows with dotted hor- 
izontal lines the frequencies of the first three unstable modes, 
i.e. S m= i, S m =3, so that S 2 can be viewed as an analog of 
the energy of a quantum mechanical particle. For each of the 
unstable modes there is therefore an inner "forbidden" region 
< r < r a , an "allowed" region r a < r < r&, a "potential 
barrier" n, < r < r c and an external "allowed" region r > r c . 
The points r a , rb and r c shown in the figure correspond to the 
£ = m = 1 mode. 

The WKB matching of the wave functions in the four re- 
gions is shown in the Appendix, where we present the deriva- 
tion of eqs. (l25T l and (|26| | below. The unstable modes are de- 
termined by requiring that 



:J VTdr= ^n+ijTr; » = 0,1,2,. 



(25) 



and 



r = 4 exp 2m 




VTdr , (26) 
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FIG. 3: Left panel: Change in the size of the ergoregion for fixed compactness fi — 0.45, several different values for 5/M (from top to 
bottom 5/M — 0.1, 0.2, 0.3 and 0.4) and increasing angular velocity Q, until the Keplerian limit Qk- A given value of Q on the vertical axis 
determines the inner and outer radii of the ergoregion, while the vertical line shows the location of the radius of the gravastar for all the models: 
the ergoregion generally starts in the interior of the gravastar and goes up to a radius exterior to the radius ri of the gravastar. Right panel: 
Same as the left panel, but for fixed thickness of the shell 5/M = 0.4 and different values for p (from top to bottom p = 0.42, 0.43, 0.45 
and 0.47). In this case we do not show the surface radius, since each model has a different radius r2. 



where 



,A-i 



T 



(27) 



Equation (l25l l is the classical Bohr-Sommerfeld rule and de- 
termines S, while eq. d26l i gives the growing time of the insta- 
bility. The limits of the integration interval r a , r& and r c have 
the physical interpretation given above for Fig. [2] and corre- 
spond, mathematically, to the turning points as given by the 
condition V+ = £ (or T — 0) and to the beginning of the 
free allowed region as given by the condition V- = S. Note 
also that the use of the absolute value for T in the first integral 
of eq. d2oT ) is due to the fact that in the interval r& < r < r c 
(i.e. inside the potential barrier) T < [cf. eq. d23lll. 

It is customary in the evaluation of the integral in eq. ( l25l ), 
to use an analytical parabolic approximation for V+ (see 
ref. lfl2lo . When this is done, the functions A, v and V- are 
taken to be constants, with their values set to \(R), v(R) and 
V-(R), where R is the radius at which V+ has its minimum. 
This greatly simplifies the calculations and is a very useful 
strategy in general. In the case of a gravastar, however, the 
potential well in V+ is typically very asymmetric around the 
minimum of the potential, so the usual analytical parabolic 
approximation is not appropriated in this case, for it would 
introduce too large errors and it has not been used in the nu- 
merical solution of eq. (f26b . 

Collected in Table U are some typical numerical values ob- 
tained for the potential V+ of the gravastar with p = 0.45 and 
5/M = 0.4 and different values of J. With these data it is 
possible to confirm that both the size of the ergoregion (which 
extends from r m ; n , the inner boundary of the ergoregion, to 
r mal , the outer boundary of the ergoregion, in the notation of 
the table) and the depth of the potential well | V+ (R) | increase 
with J, while R remains essentially unchanged. 



TABLE I: Typical values of J/M 2 , Q/Qk, r m i n (inner boundary of 
the ergoregion), r max (outer boundary of the ergoregion), R (radius 
at which V+ is minimum) and V+(R) for the gravastar with p = 
0.45 and 5/M = 0.4. 



J/M' 2 


VL/Q.K 


rndxx/M 


rmax/M 


R/M 


MV+(R) 


1.2 


0.98 


0.500 


2.411 


1.8789 


-0.2661 


1.0 


0.82 


0.591 


2.321 


1.8797 


-0.2164 


0.8 


0.65 


0.721 


2.231 


1.8806 


-0.1668 


0.6 


0.49 


0.917 


2.146 


1.8815 


-0.1172 


0.4 


0.33 


1.225 


2.071 


1.8824 


-0.0068 


0.2 


0.16 


1.683 


1.987 


1.8834 


-0.0018 






IV. 


RESULTS 







We have first checked our results against the values ob- 
tained by Cardoso et al in ref. |2[] (for £ and r) for a typical 
model with r\ = 1.8, r2 = 2.2 and M = 1, with very good 
agreement. Based on their findings for this typical model, 
which shows a very rapidly growing instability, they have con- 
cluded that the ergoregion instability would "rule out" the pos- 
sibility of rotating gravastars to exist as alternatives to Kerr 
black holes. However, drawing a general conclusion from a 
specific example can be too restrictive. In view of this, we 
have decided to reconsider the issue of whether all of the ro- 
tating gravastars are unstable to the ergoregion instability and 
if so over what timescales. As we will show in this Section, 
by considering a much larger space of parameters we reach 
conclusions which are rather different from those of ref. J2]. 

We start by showing in Fig.[3]how the size of the ergoregion 
(taken at the equatorial plane) depends on the parameters of 
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FIG. 4: Minimum angular velocity necessary for the existence of an 
ergoregion, as a function of both /i and 8. If an ergoregion is present, 
the instability will set in for a high enough value of m. In the grey 
area at the bottom, an angular velocity SI larger than the mass shed- 
ding limit £Ik would be required for an ergoregion to develop. The 
grey area at the top shows a constraint on the possible (nonrotating) 
gravastar solutions found in ref. 1 1 ill . 







'Jmin/M 2 






0.9 
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"min(Jmin) > «K 





0.1 0.2 0.3 0.4 0.5 

5/M 



FIG. 5: Minimum angular momentum necessary for the existence 
of an ergoregion, as a function of both fi and 5. Gravastars with 
parameters (fi, 5) given below a curve labeled with some value of 
Jmin/M 2 will be stable if rotating with angular momentum smaller 
or equal to J m i n . Note that the Q is not constant along the lines, and 
increases in the direction of the region labeled Sl m i n (J m i n ) > Qk- 



the gravastar and in particular on its compactness ytt, on its 
thickness 5 and on its angular velocity SI, which we take to 
range from zero to the Keplerian (mass shedding) limit Sl^ = 

(M/rl) 1 ' 2 . 

On the left panel of Fig.[3j in particular, we report the vari- 
ation of the size of the ergoregion for models with /i = 0.45 
and different values of S, with the vertical line signing the sur- 
face of the gravastar. On the right panel, on the other hand, we 
obtain essentially the same behavior, but this time for models 
with S/M = 0.4 and different values of /i. These plots are to 
be interpreted as follows: for a given value of Sl/Sl^ on the 
vertical axis, each of the curves in the plot provides the values 
of the inner and outer boundaries of the ergoregion. It is then 
easy to see that the ergoregion becomes obviously larger for 
increasing values of SI, but also for increasing values of 8 and 
li. It also should be noted that the ergoregion is mostly (but 
not exclusively) contained in the interior of the gravastar, and 
it also extends to the exterior region. Most importantly, how- 
ever, the local minima in Fig. [3] indicate the first important 
result of this investigation: not all rotating gravastars possess 
an ergoregion. Rather, for any choice of compactness and 
thickness of the gravastar, there exists a minimum angular ve- 
locity Sl m j n above which an ergoregion of finite size develops. 
In this respect, and not surprisingly, gravastars behave like ro- 
tating compact stars. It is then a trivial consequence of the 
above result that not all rotating gravastars are unstable to the 
ergoregion instability. 

This conclusion is summarized in Fig. [4] which shows the 
minimum angular velocity necessary for the existence of an 
ergoregion, as a function of both the compactness /j and thick- 
ness 5 of the gravastar. In other words, for any couple of val- 
ues of (fi, 8) outside of the shaded regions, the corresponding 



gravastar will not possess an ergoregion if spinning below the 
value of Sl mm at that point. Note also that in the grey area at 
the bottom of the figure, an angular velocity SI larger than the 
mass shedding limit CIk would be required for an ergoregion 
to develop. Similarly, the grey area at the top of the figure 
shows the constraints on the possible (nonrotating) gravastar 
solutions found in ref. Hill . It is important to remark that, as 
clearly shown in Fig. HI the minimum angular velocity Sl m ; n 
increases with decreasing compactness of the gravastars. This 
behavior sets an additional constraint on the parameters of the 
gravastars that will be subject to the ergoregion instability: 
less compact gravastars would have to rotate with SI > Six to 
form an ergoregion and are therefore also free from the insta- 
bility. 

Interestingly, there is a non-small portion of the space of 
parameters (/i, 8) where very rapidly rotating gravastars ex- 
ist, do not possess an ergoregion and are therefore stable. We 
recall, in fact, that while black holes have their angular mo- 
mentum bounded by the Kerr limit (i.e. J/M 2 < 1), stars 
(and gravastars!) are not subject to this constraint. As a result, 
as long as they are spinning below the mass-shedding limit, 
gravastar models can be built that are stable and even have 
J/M 2 > 1. This is shown in Fig.[5j which reports the min- 
imum angular momentum necessary for the existence of an 
ergoregion, as a function of both /i and 8. Gravastars with pa- 
rameters (fi, 8) given below a curve labeled with some value 
of Jmi n /M 2 will be stable if rotating with angular momentum 
smaller than J m ; n (models on the curve are marginally stable, 
i.e. with t = oo, for that value of J m i n ). This is the second 
important result of this paper: not all ultra-compact astro- 
physical objects rotating with J/M 2 ~ 1 must be black holes. 
This conclusion is thus less restrictive than the one drawn by 
Cardoso et al. in ref. J2]. 
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FIG. 6: Maximum compactness ^ max that a stable gravastar can have 
in terms of the angular momentum J/M 2 . For a given value of J, 
gravastars with /i > ^max(J) will be unstable. The maximum value 
for the compactness is obtained for 5^0. This figure should be 
compared with the corresponding Fig. 1 of ref. fill] . 



Note that the O is not constant along the J m i n — const, lines 
but, rather, it increases in the direction of the region labeled 
^min(^min) > ■ However, as the angular momentum is 
increased, the maximum compactness allowed for the stable 
models is also reduced, thus indicating a new bound on the 
compactness of stable models. We show this additional con- 
straint on the gravastar's compactness in Fig. El which reports 
the maximum allowed compactness /i max as a function of the 
angular momentum J/M 2 . [This fi gure should be compared 
with the corresponding Fig. 1 of ref. Ill] where the constraint 
on the compactness was shown as a function of the thickness 
for nonrotating gravastars.] 

Focusing on J/M 2 ~ 1, a case which is astrophysically 
very interesting l2(il l2lll . Fig. [6] shows that it is possible to 
construct stable rotating gravastars with J/M 2 ~ 1 as long 
as their compactness is less than ^ max < 0.43. Such a com- 
pactness is clearly smaller than that of a black hole with the 
same angular momentum (i.e. /.t = 1), but much larger than 
the typical compactness for compact stars and neutron stars, 
i.e. 0.15 < [i < 0.2. (We recall that the maximum compact- 
ness for for perfect fluid nonrotating spheres is given by the 
Buchdahl-Bondi limit fx = 4/9 ~ 0.44 JUIH). It is still un- 
clear whether astronomical electromagnetic observations will 
(ever) be able to distinguish a stable rapidly rotating gravastar 
with [i ~ 0.43 from a rotating black hole with J/M 2 ~ 1 
(see the discussion in ref. 12411 ). Yet, the constraints emerging 
from Fig.|6]provide other means (besides the measurement of 
quasi-normal modes Hill ) in which astronomical observations 
could be used to distinguish (rotating) gravastars from (rotat- 
ing) black holes. 

So far all of the considerations made were on the general 
properties of rotating gravastars. However, to fix the ideas 
and also provide some reference numbers on the frequencies 
and timescales for the ergoregion, we now discuss in more 



detail one typical case which we will assume to be our ref- 
erence model (cf. Table U]). More specifically, we consider a 
gravastar with mass M = IOMq and inner and outer radii 
n = 1.8M ~ 27 km, r 2 = 2.2M ~ 32 km (or, alternatively, 
fi ~ 0.45, 5/M = 0.4). Its mass-shedding spin frequency 
will be vjc = Hk/2tt ~ 990 Hz. Such a gravastar will be 
stable (because without ergoregion), for J/M 2 < 0.13 (or 
£1/£Ik < 0.10). Conversely, if the gravastar is set to rotate at 
a higher rate, namely with J/M 2 ~ 0.22 (or Q/Q K ~ 0.18), 
an ergoregion will be present and the instability will develop 
over a timescale t ~ 10 17 s, with I = m = 4 being the lowest 
unstable mode. Note that such a timescale is comparable with 
the Hubble timescale and thus the rotating gravastar will be 
physically stable although mathematically unstable (a similar 
result holds true also for compact uniform-density stars lfl2"tl ). 
In practice, it is necessary to spin the same gravastar up to 
J/M 2 ~ 0.61(orfi/O A ' ^ 0.50) in order for the instability to 
develop on a timescale of the order of 1 s, with the £ = m = 1 
being the lowest unstable mode. 

A more complete picture of the real and imaginary parts 
of the eigenfrequencies for our representative rotating gravas- 
tars is shown in Fig. [7] where we have extensively explored 
the parameter space 5, J). More specifically, the top 
row shows in the left panel the timescale r of the instabil- 
ity (with I = m = 1) as a function of the angular mo- 
mentum J of the gravastar, for fixed thickness of the shell 
5/M = 0.4 and different compactnesses [i (from left to right 
// = 0.47, 0.45, 0.44 and 0.43). As a reference the horizontal 
line shows a timescale of Is for gravastars with M = IOMq. 
Similarly, the right panel of the top row gives the frequency 
£ as a function of J for the same cases in the left panel. 
The bottom row, on the other hand, shows the same quanti- 
ties as the top one, but for gravastars with fixed compactness 
fi = 0.45 and varying thickness of the shell 5 (from left to 
right 5/M = 0.6, 0.5, 0.4, 0.3, and 0.2). 

Note that the timescale for the growth of the instability in- 
creases over-exponentially with decreasing angular momen- 
tum J, and this is why a small spin down of the gravastar is 
sufficient to make it "effectively stable" even if it possesses an 
ergoregion and is mathematically unstable. Note also that the 
frequency £ increases only linearly with decreasing J. Fi- 
nally, to aid those interested in reproducing our results and 
for code-testing purposes, we have collected in Table HT1 some 
numerical values for the first unstable mode relative to our ref- 
erence gravastar (cf. Tabled. Clearly, the values of r are very 
sensitive to J, but the big difference in r seen between the 
gravastar with J/M 2 = 0.4 and the one with J/M 2 = 0.2 
is also due to the fact that in the first case the lowest unstable 
mode is still I = m = 1, while in the second case it is already 
I = m = 5. 

The concluding remark of this Section should be one of 
caution. All of our treatment is based on the slow-rotation 
approximation, yet we have stretched it to compute also mod- 
els with 0,/VIk = 1. Furthermore, while the WKB formulas 
we used were derived in the high-m limit, we have effectively 
used them also for very low values of m and found, in par- 
ticular, that most rotating gravastars are unstable already for 
I = m = 1 mode, as a result of their high compactness. Al- 
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FIG. 7: 7bp row, left panel: Timescale r of the instability (with I — m = 1) as a function of the angular momentum J of the gravastar, for 
fixed thickness of the shell 5/M = 0.4 and different compactnesses /i (from left to right /1 = 0.47, 0.45, 0.44 and 0.43). The horizontal line 
shows a timescale of Is for gravastars with M = IOMq. Top row, right panel: Frequency E of the instability as a function of J for the same 
cases in the left panel. Bottom row: The same as in the top row, but for gravastars with fixed compactness /j, = 0.45 and varying thickness of 
the shell 6 (from left to right 5/M = 0.6, 0.5, 0.4, 0.3 and 0.2). 



though our approach is not dissimilar to the one made in re- 
lated works J3[I3l, it is important to underline that we expect 
our estimates to be accurate only for slowly rotating gravas- 
tars. On the other hand, we also believe that the qualitative 
(and possibly quantitative) picture derived here will remain 
unchanged also when a more sophisticated analysis is per- 
formed. 



V. CONCLUSIONS 

Motivated by recent work on this subject [Q], we have in- 
vestigated the ergoregion instability in rotating gravastars, ex- 
ploring a large space of parameters and taking into account 
the limits on the thickness of the matter shell and on the com- 
pactness. While we confirm the results of Cardoso et al. yfl 
for the models they have considered, we also draw two con- 
clusions which are less restrictive than theirs. Firstly, we find 
that models of rotating gravastars without an ergoregion (and 
therefore stable) can be constructed even for extreme rotation 



TABLE II: Typical values of J/M 2 , r a , rt, r c , E and r for for the 
first unstable mode of the gravastar with [i = 0.45 and S/M = 0.4. 
Note that all values but the last one refer to the I = m = 1 mode. 

J/M 2 r a /M r b /M r c /M MS r/M unstable mode 

1.2 1.24 2.048 4.31 -0.20 1.24 x 10 3 i = m = l 

1.0 1.24 2.052 5.76 -0.15 2.35 x 10 3 i = m = l 

0.8 1.24 2.056 8.88 -0.10 5.96 x 10 3 t = m = l 

0.6 1.25 2.061 18.1 -0.052 2.73 x 10 4 £ = m = l 

0.4 1.25 2.066 350 -0.0027 1.12 x 10 7 i = m = \ 



0.2 



1.70 1.982 592 -0.0017 1.19 x 10 30 i = m 



rates, namely for models with J/M 2 ~ 1. Hence, not all 
rotating gravastars possess an ergoregion. Secondly, because 
stable gravastar models with J/M 2 ~ 1 can be constructed, 
we conclude that not all ultra-compact astrophysical objects 
rotating with J/M 2 ~ 1 must be black holes. 

Besides clarifying these two important aspects of rotating 
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gravastars, our analysis also helps to further constrain the 
properties of these ultra-compact objects. Building on our ini- 
tial work ITTll . in fact, we have computed an additional con- 
straint on the maximum compactness of a gravastar which is to 
be stable to the ergoregion instability. Such a maximum com- 
pactness is still much larger than that of typical neutron stars 
but also smaller than that of black holes with J/M 2 ~ 1. This 
should help in the important effort of distinguishing (rotating) 
gravastars from (rotating) black holes. 
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APPENDIX 

We review the derivation of the WKB formulas dZSb and 
used in this paper, following closely the treatment given 
in ref. B12I1 . We have chosen to include this Appendix in the 
paper, even though the usual WKB approximation is standard 
in quantum mechanics textbooks, because the case in ques- 
tion is somewhat more complicated than the usual textbook 
exercises. The equation to be solved is [cf. eq. d22b l 



ip irr + m 2 T(r, £)V> = , 



(A28) 



which has four different regions with distinct physical behav- 
ior, as described in Section[ni] We are interested in finding the 
purely outgoing modes of this equation. These are given as 
poles of the scattering amplitude S = C ou t/C- m , where C ou t 
and Cin are functions of the complex frequency a (S = a jm 
and we assume a e lcrt dependence for ip) and denote the am- 
plitude of the outgoing and incoming waves at infinity, respec- 
tively. In a pole of S, we will have C m = and C ou t ^ 0. For 
the modes with small imaginary part, 5R(cr) ^ 3(c) it will be 
sufficient to determine the eigenvalue a approximately as a 
real frequency (on the real axis). In this spirit we define an 
auxiliary function S(a) = [S(a*)}*, such that if S has a pole 
at <7p, then S has a pole at tr*, where (*) denotes complex 
conjugation. 

If we restrict ourselves now to a e R, we will have conser- 
vation of energy and 

\S\ = 1 => S(a)[S(a)r = 1 => S(a) = [5(a)]- 1 , 

(A29) 

where the last result is obtained because a is real. As the 
relation above is valid for all real values of a, it is valid every- 
where S and S are analytic functions. Therefore, returning to 
the complex plane, we can state that if S has a pole at cr*, then 
S has a zero at a* (and still a pole at a p , which we assume to 



be simple), and can be written approximately as 

S = e 2i So ZZfl g = e 2i ^-^-i/T (A3()) 

a — Up a — (T r + l/r 

Now we turn to the specific form of the solutions of 
eq. (|A28t in the four different regions as discussed in Sect. [HI] 
region I or inner "forbidden" region, region II or "allowed" 
region, region III or "potential barrier", region IV: r > r c 
or external "allowed" region {cf. Fig. [2j. The connection for- 
mulas for the wave function to the left and right of a turning 
point can be found in many standard textbooks on quantum 
mechanics {e.g. JHI) and amount to 
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for a turning point x 
and 



a to the right of the classical region, 



1 



: CXp 



: exp 




kdx 



kdx 



(A33) 



J)(A34) 



for a turning point x = b to the left of the classical re- 
gion. In this notation, the equation to be solved is d 2 ip/dx 2 + 
k 2 (x)ip — or d 2 ip/dx 2 ~ n 2 (x)ip = 0, with k 2 , k 2 > 0. We 
will need this formulas in what follows, to connect the wave 
function in the different regions. The form of the radial func- 
tion ip( r ) m region I is determined by the regularity condition 
at the center {ip must vanish at r = 0), 



01 



„l/2| T |l/4 



exp — m 



T\dr 



(A35) 



Using now the connection formula ( |A33t , we can obtain the 
form of the radial wave function in region II (connecting 
through the first turning point r = r a ) 



4>n 



+ r l/2 T l/4 



exp lm 



Tdr 



exp —l m 



Tdr 



(A36) 



where £ 
solution in regions III and IV as 

C 2 



7r/4. We can now write the 



in 



r l/2| T |l/4 

c 3 



exp 



exp 
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r l/2 T l/4 



exp 




(A37) 



(A3 8) 
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Combining the expressions for tpm and ipiy with the con- 
nection formulas dA3 lb -( |A34t , we can relate the amplitudes 
of the waves before and after crossing the potential barrier. 
This is done as in a standard scattering exercise, demand- 
ing the continuity of the wave function and its first derivative 
across the turning points r?, and r c . The result, as stated in 
ref. i25ll . is remarkably simple and best expressed in matrix 
notation: 




2r, 



j_ 



2rj 




C 4 \ 

c 5 J 

(A39) 

We can now identify the 



where 77 = exp (m \/\T\ dr 

amplitudes C4 and C5 with the amplitudes C m and C ou t of the 
incoming and outgoing waves at infinity that we have intro- 
duced in the definition of S. If E is negative, the case in which 
we have unstable modes, then C4 = C ou t and C5 = C- m . In- 
verting now eq. JA39t in order to obtain C4 and C5 in terms 
of C\, they can be substituted in the definition of S to give 



Q _ (4?? 2 + l)e^ - i( V - l)e~* 
~ ~Cl ~ i(V - l)e* + (V + l)e-'< 



(A40) 



As we have discussed above in the beginning of this Ap- 
pendix, we are interested in finding the purely outgoing 
modes, which are poles of the scattering amplitude S. Also, 
we made the assumption that this poles will occur for complex 
frequencies er p which lie on the complex plane, but very close 
to the real axis. But the imaginary part of <r p will only be 
very small if the "barrier penetration" integral 77 is very large. 
Taking the limit of S as 77 — > 00, we have 



-K 



-i for 77 — > 00, 



(A41) 



unless we have e 1 ^ — ie ^ = 0, in which case we can write S 
as 



4 7? 2( ie iC + e -iC)_( ie iC_ e - i C) ( ' 



Therefore we can see that S will have a resonance at a fre- 
quency <r n near the frequency for which e 1 ^ — ie -1 ^ = and 
thus £ = wr + 7r /4, with n an integer. We can now write £ as a 
series expansion around cr n , £(er) ~ n n + tt/4 + a n (<7 — <x„), 
where 



da 



Tdr 



(A43) 



Substituting now £(er) into eq. dA40| ), we obtain at first order 
in (a — <j n ) and dropping the index n on a n 



S 



-a(a - a n ) + I/477 2 + i [a{cr - a n ) + 1/47/ 2 ] 
-a(a - a n ) + I/477 2 - i [a{a - a n ) + I/477 2 ] ' 

(A44) 



which can then be rewritten in the form of eq. ( |A30t , 



S = 



.a — a n — i/Arfa 



(A45) 



Combining the definition of C, and its expansion and taking 
cr = a n we finally obtain eq. d25l l and comparing eqs. JA30b 
and dA45t , we obtain eq. ( f26T >, both presented in SectionllTII 
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